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Lagrangian dispersion and heat transport in convective turbulence 
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Lagrangian studies of the local temperature mixing and heat transport in turbulent Rayleigh- 
Benard convection are presented, based on three-dimensional direct numerical simulations. Contrary 
to vertical pair distances, the temporal growth of lateral pair distances agrees with the Richardson 
law, but yields a smaller Richardson constant due to correlated pair motion in plumes. Our results 
thus imply that Richardson dispersion is also found in anisotropic turbulence. We find that extremely 
large vertical accelerations appear less frequently than lateral ones and are not connected with rising 
or falling thermal plumes. The height-dependent joint Lagrangian statistics of vertical acceleration 
and local heat transfer allows us to identify a zone which is dominated by thermal plume mixing. 

PACS numbers: 47.55.pb,47.27.tc 
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Introduction. The interplay of turbulent temperature 
and velocity fields in thermal convection is present in 
many systems reaching from astro- and geophysical flows 
0, Q to chemical engineering and indoor ventilation Q . 
Almost all experimental and numerical studies on ther- 
mal convection have been carried out in the Eulerian 
frame of reference where the velocity and the tempera- 
ture are analysed at fixed spatial positions and the mean 
turbulent heat transfer is measured with increasing pre- 
cision for ever larger Rayleigh numbers [1, H, @. Not 
much is however known about the local mechanisms of 
the transport of heat and momentum. Heat transport 
fluctuates locallly since it is arranged in coherent ther- 
mal plumes. They detach permanently from the thermal 
boundary layers, carry blobs of hotter or colder fluid into 
the bulk and cause an efficient temperature mixing. The 
complementary Lagrangian description in which turbu- 
lent fields are monitored along the tracks of fluid parcels 
0, H |, M, ED, 03, El could thus provide exactly 
the missing insights on the local mixing and transport 
mechanisms in convective turbulence. For example, the 
frequency and seasonal variability of thermal plumes af- 
fects the transport and concentration of phytoplankton 
in the upper ocean and has an impact on the carbon cycle 
[l5L . Vertically transported graupel particles or wa- 
ter droplets in deep-convective cumulus clouds determine 
how much precipitation falls out [IT], [l8| . 

In this Letter, we study the Lagrangian mixing of tur- 
bulent Rayleigh-Benard convection in three-dimensional 
direct numerical simulations. Two aspects are in the fo- 
cus of the present investigation. Firstly, we provide a 
detailed analysis of the acceleration and dispersion prop- 
erties of Lagrangian tracers in convection. Caused by the 
interplay of buoyancy and gravity, Rayleigh-Benard tur- 
bulence is inhomogeneous and anisotropic. Consequently, 
significant differences between the vertical and lateral 
transport properties are found. Secondly, we will take 
a "Lagrangian fingerprint" of thermal plumes. We there- 
fore explore the statistics of local heat transfer and its 
relation to the local acceleration along the Lagrangian 
tracer tracks, all this as a function of height. This al- 




FIG. 1: (Color online) Instantaneous snapshot of two isosur- 
faces of the total temperature T (top=cold, bottom=hot). 
Thermal plumes which detach from both planes can be ob- 
served. Streamlines of the velocity field are also shown. 



lows us to identify the regions in the convection cell 
which are dominated by the plumes and in which conse- 
quently enhanced vertical temperature mixing is present 

. Our studies extend recent Eulerian investigations on 
the plume structure [n| [2(| and an experiment with neu- 
trally buoyant smart temperature sensors which was able 
to probe Lagrangian bulk properties in high-Rayleigh- 
number convection on scales down to 17 times the ther- 
mal boundary layer thickness [2lj . 

Numerical simulations. The Boussinesq equations for 
an incompressible flow and the advection-diffusion equa- 
tion for the temperature field, are solved by a standard 
pseudospectral method for the three-dimensional case. 
The equations are given by 



V u = 0, 



(1) 



du 



^ + (u-V)u = -Vp + vW 2 u + ag6e z , (2) 
at 



?l + {u-V)6 = nV 2 6 + u z ^-. 



(3) 



Here, u is the turbulent velocity field, p the pressure 
field and 9 the temperature field. The system parameters 
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FIG. 3: (Color online) Probability density function (PDF) of 
the acceleration vector components. Components are normal- 
ized by their corresponding root-mean-square values (rms). 
The skewness of all three distributions vaires between -0.2 
and -0.09. The dashed vertical lines mark the characteristic 



acceleration u z 



(v) 



FIG. 2: (Color online) Particle pair dispersion in turbulent 
convection, (a): Pair dispersion ([A(£) — A(0)] 2 )z, as a func- 
tion of time t. The vector A is the distance vector between 
both tracer positions within a pair. The vertical, lateral and 
total dispersion are compared. The height of the convection 
cell, the Kolmogorov time, t v , the lateral Lagrangian integral 
time, Tl, are also indicated, (b): Probability density func- 
tion (PDF) of the lateral dispersion taken at four times, see 
symbols in (a). For comparison, we include the Gaussian pre- 
diction by Batchelor [23] and the Richardson prediction [22l ] 
(with different fit factors). The inset shows the Lagrangian 
velocity autocorrelation functions for the lateral, Ci a t{r) and 
vertical , C ve r(j) components. The vertical dashed line cor- 
responds with the travel time H/u z ,rms- 

are: gravity acceleration g, kinematic viscosity v, thermal 
diffusivity k, vertical temperature gradient AT/H, and 
thermal expansion coefficient a. The total temperature 
field T is decomposed into a linear background profile 
and a fluctuating field 

AT 

T(x,t) = -—(z-H/2)+6(x,t). (4) 

Since T = ±AT/2 at boundaries z — and z = 
H(= 7r), the boundary condition 8 = follows. Here, 
AT > 0. The dimensionless control parameters are the 
Prandtl number Pr = v / n = 0.7, the Rayleigh number 
Ra = agH 3 AT/(vn) = 1.2 x 10 8 , and the aspect ratio 
r = L/H = 2 with L = L x = L y = 2tt. Figure 1 il- 
lustrates a typical snapshot of the flow and temperature 
fields. In the vertical direction z, free-slip boundary con- 
ditions are used. Laterally in x and y, periodic boundary 
conditions are applied. The simulation box is resolved 
by 512 x 512 x 257 grid points which results in a spectral 



resolution of k max r)K = 2.2. The Taylor microscale The 
Reynolds number is R\ = 135 and the Nusselt number 
Nu = 58.05 ± 0.98. Time-stepping is done by a second- 
order predictor-corrector scheme. We track N = 10 6 La- 
grangian tracers simultaneously with the flow. Initially, 
they are arranged in pairs with a separation of 1.5 grid 
cell widths. Intergrid velocities are obtained by trilinear 
interpolation. 

Particle dispersion. Figure 2a reports our findings on 
the pair dispersion which is given by R 2 (t) = ([A(t) — 
A(0)] 2 )l. The statistics is taken over N/2 tracer pairs 
and denoted as Each pair is described by the par- 

ticle distance vector A = X2 — x\. Interestingly, af- 
ter the initial ballistic growth with t 2 a crossover to a 
Richardson- like dispersion with R 2 (t) = C2(e)i 3 is ob- 
served, similar to isotropic fluid turbulence @, [ll|, |22| . 
Here, (e) is the mean energy dissipation rate. In con- 
trast to the isotropic case, the Richardson constant is 
however reduced by two third and is given by C2 = 0.16 
in the present study 0, H3 ■ Recall, that the turbulence 
is inhomogeneous with respect to the vertical direction. 
Therefore, we decompose the distance vector into a lat- 
eral and vertical part, A = r xy + r z e z . While the vertical 
distances cannot grow beyond the cell height H, the lat- 
eral one takes over and makes up the whole large-time 
growth of R 2 (t) including the Richardson scaling. Even- 
tually, the lateral pair motion becomes decorrelated and 
changes to the diffusive limit with R 2 (t) ~ t (see Fig. 
2a). 

In Fig. 2a, the viscous Kolmogorov time r, ; = y/v/{e) 
and the lateral Lagrangian integral time scale Tl are 
given. Both are separated by less than an order of mag- 
nitude for the Reynolds number accessible here. Fur- 
thermore, the i 3 -regime extends beyond Tl- Figure 2b 



3 



shows the probability density function (PDF) of the lat- 
eral distance at five different times. The PDFs collapse 
to a Richardson-like shape and clearly deviate from the 
Gaussian limit for times around T^. They cross over to 
Gaussian shape for larger times. A Gaussian distribu- 
tion, as predicted by Batchelor for the Lagrangian iner- 
tial range, is not supported by the present study [23| . 

The inset of Fig. 2b unravels another specifics in con- 
vection. The Lagrangian velocity autocorrelation func- 
tions of the vertical and lateral components, C ver (r) and 
Ci a t(r), which are calculated along the tracer tracks are 
shown. The strong anticorrelation of the vertical part 
is due to the constraints in vertical motion, e.g. when 
the tracers (and the plumes) hit the top and bottom 
planes. The pronounced minimum in the vertical cor- 
relation is close to a characteristic vertical travel time 
~ H/u z , rms , where u z , rms — (u 2 ,) 1 / 2 . The integration 
of C ver (r) with respect to time results in a vertical La- 
grangian integral time which is ps Tl/6. While the 
single particle motion decorrelates on a shorter time scale 
as in isotropic turbulence, the lateral pair motion seems 
to be stronger correlated. The latter is due to tracer 
pair motion in plumes. This circumstence rationalizes 
the smaller Richardson constant C2 and why R 2 (t) ~ i 3 
seems to continue for times t > Tl. Our studies demon- 
strate, that Richardson dispersion can also be found in 
anisotropic turbulence. They could thus shed light on 
why Richardson-like dispersion is frequently observed, 
but the constant C2 strongly scatters around the ana- 
lytically predicted value. 

Lagrangian plume detection. The results on the par- 
ticle pair dispersion suggest to take a closer look to the 
vertical accelerations and how their amplitudes are cor- 
related with the thermal plumes. In Fig. 3, we show 
first the PDFs of all three acceleration components. As 
expected, the PDFs of both lateral components, a x and 
a y , collapse onto each other. Surprisingly, they yield fat- 
ter tails compared to the remaining vertical acceleration 
component. The characteristic amplitude for the latter 
can be estimated by \a z \ ~ u z ^ rrns /T^ which is indi- 
cated in Fig. 3. The largest vertical acceleration ampli- 
tudes result then to about 10 times this estimate. The 
overall symmetric and stretched exponential shape re- 
veals the strong intermittency of Lagrangian motion as 
in isotropic fluid turbulence [T2, EH ■ The differences in 
lateral and vertical accelerations suggest that the rise of 
buoyant plumes is less gradual than the lateral vortical 
motion. 

It is known that rising or falling plumes can be iden- 
tified by a positive product of vertical velocity and 
temperature fluctuation, u z T' > where T'(x,t) = 
T(x,t) — (T)a(z). Here, (-)a denotes ensemble averages 
at fixed height z. The product u z T' is directly related to 
the (Lagrangian) Nusselt number [2l[ which is given by 
Nul = 1 + H/(kAT)u z T'. Figure 4c shows how prob- 
able the positive and negative amplitudes of u z T' as a 
function of height z are. We see, that the distribution 
is always asymmetric with respect to the axis u z T' = 0. 
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FIG. 4: (Color online) Lagrangian heat flux as a function 
of height, (a) Vertical slice cut through the instantaneous 
temperature field T. (b) Corresponding product field u z T' ' . 
(c) Lagrangian probability density function (PDF) of u z T' 
as a function of vertical tracer position z. The height z = 
H/3 corresponds with 38 times the thermal boundary layer 
thickness St- The regions between both thick horizontal lines 
at H/3 and H/16 as well as 2H/3 and 15if/16 are considered 
as the mixing zones. Color legend is here in decadic logarithm. 



This reflects the net transfer of heat from the bottom 
to the top plane. The broad support of the PDF is 
due to the fluctuations of the local heat transfer. The 
joint PDF shows two features. Firstly, the support of 
the PDF for u z T' > is broadest at about 18 times 
the thermal boundary layer thickness St (i.e. z sw if/6). 
Secondly, the support of the PDF for u z T' < peaks at 
if/16 and decreases to about if/3. In order to shed more 
light on the structural features which are associated with 
this statisitics, we show in Fig. 4 the temperature T (a) 
and the product u z T' (b). Clearly, thermal plumes are 
found were u z T' > 0. However, one observes also larger 
u z T' < in the vicinity of plumes. Due to the incom- 
pressibility of the flow, strong upwellings are accompa- 
nied by neighboring downwellings which reflects in the 
statistical distribution in the bottom panel of Fig. 4. By 
combining both aspects, we conclude that the heat trans- 
port and mixing are dominated by thermal plumes in the 
zone between if/ 16 % z < if/3. Such plume mixing re- 
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FIG. 5: (Color online) Joint Lagrangian PDF of vertical ac- 
celeration component a z /a Zirms and the product u z T' . Color 
coding is again in decadic logarithmic units as indicated by 
the color bar. The regions for rising hot plumes and falling 
cold plumes are indicated by arrows. 

gion has been suggested by Castain g et al. [H and was 
determined in an Eulerian analysis [24j . Compared to 
[24j , the present mixing zone is smaller in extension, but 
again significantly larger as the thermal boundary layer 
thickness. This explains why the smart sensor experi- 
ment by Gasteuil et al. [2j} could detect successfully the 
large fluctuations of the heat transfer despite the large 
probe diameter. 

Finally, the combination of the results on accelera- 
tion and local heat transfer result in the joint statistics 
which is shown in Fig. [5] In order to highlight the 
statistical dependence of both fields, we normalize the 
joint PDF by the single quantity PDFs, U{a z ,u z T) = 



P(a z ,u z T')/[P(a z )P(u z T')}. The pronounced asymme- 
try with respect to the axis u z T' — reflects again the 
net transfer of heat, similar to Fig. 4. It can be seen that 
larger accelerations are associated with rising and falling 
plumes as well as with resulting down- and upwellings, 
respectively. The latter are found in the corresponding 
diagonal quadrant of the plane. Maxima of II(a z , u z T') 
are detected in the outer regions of the support, but not 
for the strongest acceleration. The figure illustrates thus 
also which fraction of the tail of the acceleration PDF is 
most probably associated with the plume motion. 

Summary. Due to vertical buoyancy driving of convec- 
tive turbulence, Lagrangian tracer motion is found to be 
strongly anisotropic in contrast to classical fluid turbu- 
lence 0, [ll|. On the one hand, the vertical Lagrangian 
transport is found significantly shorter-in-time correlated 
than the lateral one. On the other hand, rising thermal 
plumes seem to enforce a more coherent lateral motion 
which results in a smaller Richardson constant for the lat- 
eral Richardson pair dispersion which has been detected 
here in an anisotropic flow. Vertical accelerations are di- 
minished compared to lateral ones and their largest am- 
plitudes are not connected with rising and falling plumes. 
Based on the Lagrangian monitoring of heat transfer, the 
plume mixing zone has been identified complementary to 
previous Eulerian analysis. It is believed, that the pre- 
sented features can be carried over qualitatively to the 
case of no-slip boundaries. This particular investigation 
and trends with increasing Rayleigh number are part of 
our future work. 
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